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Abstract 

We perform a numerical investigation of the fields evolution in the super- 
symmetric inflationary model based on radiative corrections. Supergravity 
corrections are also included. We find that, out of all the examined initial 
data, only about 10% give an adequate amount of inflation and can be con- 
sidered as "natural". Moreover, these successful initial conditions appear 
scattered and more or less isolated. 
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A couple of years ago Linde has proposed, in the context of non-supersymmetric 
(SUSY) grand unified theories (GUTs), a clever inflationary scenario based on a coupled 
system of two real scalar fields one of which is a gauge non-singlet. This field remains 
displaced from its vacuum value during inflation and provides the vacuum energy which 
drives inflation while the other, which is a gauge singlet, is the slowly varying field during 
inflation. The main advantage of this hybrid inflationary scenario is that it can reproduce, 
in contrast to previous realizations of inflation, the observed temperature fluctuations in the 
cosmic background radiation (CBR) with natural values of the coupling constants. However, 
inflation terminates abruptly and is followed by a "waterfall" regime during which topological 
defects are easily produced if they are predicted by the symmetry breaking associated with 
inflation. 

Recently two variants of the above scenario, both sharing its naturalness, were 

proposed in the context of SUSY GUTs. One of them || relies on the lowest order non- 
renormalizable contribution to the superpotential. This scenario has been called smooth 
hybrid inflation since the termination of the inflationary period is not as abrupt as in Linde's 
scheme. Rather, the entrance of the system into the oscillatory phase is quite smooth. 
Moreover, the system, already from the beginning of inflation, follows a particular valley 
of minima that leads to a particular point of the vacuum manifold and, thus, there is no 
production of topological defects in this case. This is certainly an advantage of this scenario 
since it makes it applicable even to GUT models where the symmetry breaking associated 
with inflation does imply the existence of (possibly) cosmologically disastrous topological 
defects like domain walls or magnetic monopoles. The main advantage though of smooth 
hybrid inflation is that it can be associated to the SUSY GUT scale Msg — 2 x 10 16 GeV 
consistent with the unification of the minimal supersymmetric standard model (MSSM) 
gauge couplings. A possible shortcoming of this scenario is that the relevant part of inflation 
takes place at relatively high values [(1 — 3) x 10 17 GeV] of the inflaton field and, thus, 
supergravity and superstring corrections may not be so easily tamable. 

The second SUSY variant fU of the hybrid inflationary scenario utilizes the one loop 
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radiative corrections to the inflationary potential. These corrections are non-trivial since 
SUSY is broken at the inflationary trajectory in field space. The termination of inflation is 
as in the non-SUSY Linde's scenario and topological defects can be copiously produced. It 
is then important to make sure that no cosmologically dangerous topological defects come 
into existence as a result of the symmetry breaking associated to inflation. Moreover, as it 
turns out, the scale at which this symmetry breaking takes place is about five times smaller 
than Msg ■ This can be considered as a disadvantage of the scheme since it forces us to 
abandon the successful unification of the MSSM gauge coupling constants. In this case, 
inflation may be chosen to take place during an intermediate symmetry breaking with no 
magnetic monopoles or domain walls associated with it. The important advantage of this 
scenario is that the relevant part of inflation takes place at values of the inflaton field of 
order 10 16 GeV which are low enough for supergravity and superstring corrections to be 
under relatively good control ||]. 

For an inflationary scenario to be considered fully successful, one has to show that it is 
obtainable for a wide class of " natural" initial values of the fields and their time derivatives. 
In ref. |jj this issue was investigated for smooth hybrid inflation. The result was striking 
and unexpected. It was found that, for almost all the examined initial conditions (except for 
a narrow transition region), we get inflation with an adequate number of e-foldings. This is 
certainly an important advantage of this inflationary scenario. The purpose of this work is 
to discuss the question of initial conditions for the second SUSY hybrid inflationary scenario 
which is based on radiative corrections. In particular, we will try to see if there exists a wide 
class of "natural" initial conditions, i.e., a class of comparable initial values of the fields with 
vanishing time derivatives for which the system falls at the bottom of the valley of minima 
so that its subsequent evolution along this valley produces an adequate amount of inflation. 
To this end, we solve numerically the evolution equations of the system for various sets of 
initial conditions. Our analysis also includes cases with all initial field values being much 
smaller than the Planck scale, where the results are expected to be less affected by replacing 
global by local SUSY. At the end, non-vanishing initial time derivatives of the fields as well 
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as supergravity corrections were also taken into account. 

The SUSY inflationary scenario based on radiative corrections can be realized in the 
context of a SUSY GUT model with a (semisimple) gauge group G of rank > 5 which 
breaks spontaneously directly to the standard model (SM) group Gs at a scale M. The 
symmetry breaking of G to Gs is obtained through a superpotential which includes the 
terms |§ 

W = kS(-M 2 + 00). (1) 

Here, 0, is a conjugate pair of left-handed SM singlet superfields which belong to non-trivial 
representations of the gauge group G and reduce its rank through their vacuum expectation 
values (vevs) and S is a gauge singlet left-handed superfield. The superpotential in eq.(|l|) is 
in fact the only renormalizable superpotential consistent with a continuous U(l) R-symmetry 
under which W — > e %e W , S — > e S, 00 — > 00 ■ The potential obtained from W in eq. ([!]), 
in the SUSY limit, is 

V = k 2 | M 2 -00 | 2 +k 2 | S | 2 (| | 2 + | | 2 ) + D -terms ■ (2) 

Vanishing of the D-terms is achieved along the D-flat directions where |0| = \<p\ . The SUSY 



vacuum 



< S >= 0, < (f) >< (f) >= M 2 , | < (f) > | = | < (f) > | (3) 

lies on the particular D-flat direction 0* = 0. Restricting ourselves to this direction and 
performing appropriate gauge and R-transformations we can bring the complex S, 0,0 fields 
on the real axis, i.e., S = = 0= |x, where a and x are rea l scalar fields. The 
potential in eq.(|2]) then takes the form 

V( X ,a) = K 2 (M 2 -^) 2 + ^^ ■ (4) 

The vacua lie at < x >— ±2M, < o >= 0. We observe that the potential possesses an 
exact flat direction at x = with V(x = 0, o") = k, 2 M 4 . The mass squared of the field x 



along this flat direction is given by m 2 = —k 2 M 2 + \k 2 <j 2 and remains non- negative for 
a > a c = v2M. This implies that, at x = an d ° > a c, we obtain a flat-bottomed valley 
of minima which, as it stands, is not suitable for inflation. However, due to SUSY breaking 
along the bottom of this valley resulting from the non-zero value of V(x = 0, a), there exist 
non-trivial radiative corrections to the effective potential so that this valley acquires the 
negative slope which is necessary for inflation. The one loop corrected effective potential 
(along the direction a > o c , x = ) is given by |§ 



Veffia) = kV + _^ [2M * h, + (i S - ^) In (l - a^) + 



l , A 2 / 2M 2 

-ex + M | In 1 + 



2" ' """ J ~ \~ ' a 2 
If cr 3> cr c , V^//(cr) in eq.(|J) reduces to the simpler form 



(5) 



V eff (a > a c ) w k 2 M a 
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A region of the universe, where x an d cr happen to be almost uniform with negligible 
kinetic energies and with values close to the bottom of the valley of minima, follows this 
valley in its subsequent evolution and undergoes inflation. The temperature fluctuations of 
CBR produced during this inflation can be estimated to be [|J 

5T /iVq f M x 2 



T V 45 \M ( ' ^ 

where M c = M p /V8tt, M p = 1.22 x 10 19 GeV is the Planck mass, and N Q ~ 60 is the 
number of e-foldings of our present horizon scale during inflation. The Cosmic Background 
Explorer (COBE) result, ST/T ~ 6.6 x 10~ 6 , can then be reproduced with M = 5.82 x 10 15 
GeV. Using eq.(|6|) one finds 

where a Q is the value of a when the scale which evolved to our present horizon size crossed 
outside the de Sitter horizon during inflation. Inflation continues until a = a c where 
it terminates abruptly and is followed by a "waterfall", i.e., a sudden entrance into an 



oscillatory phase about a global minimum. This is checked by noting the magnitude of 
the quantities e = (V /V) M c 2 /2, 77 = M 2 V"/V, where the prime refers to derivatives with 
respect to o . Since the system can fall into either of the two available global minima with 
equal probability, topological defects can be easily produced if they are predicted by the 
relevant symmetry breaking of the particular particle physics model one is employing. 

We will now try to specify the initial conditions for the a and x fields which lead to the 
above described inflationary scenario. In other words, we will try to identify initial conditions 
for which the system falls at the bottom of the valley of minima of the effective potential 
so that its subsequent evolution along this valley produces an adequate amount of inflation. 
We assume that, after " compactification" at some initial cosmic time, a region emerges in 
the universe where the scalar fields a and x an d their time derivatives happen to be almost 
uniform and the energy density < M£. (The initial values of a and x together with their 
time derivatives can always be transformed by appropriate gauge and R-transformations to 
become positive). The radiative corrections in eqs.(|5]) and (|6|), though crucial for inflation, 
play no essential role in the issue of initial conditions and, thus, they will be ignored. In 
order to keep the formalism as clear as possible, we shall from now on use dimensionless 
variables for all relevant quantities as follows: 

Y = , cr = , t = KM c t, M = , H = , = — — -, V y, a) = — - — — , 9 

A M c ' M c ' M c ' kM c ' y k 2 M^ VA ' ; k 2 M c 4 ' v ; 

where H is the Hubble parameter and g the energy density. The evolution of the system in 
this region is governed by the following equations of motion 

i + 3H X - x(M 2 -^) + ^a 2 X = 0, (10) 

'a + 3Ha + ^x 2 a = , (11) 

where overdots denote derivatives with respect to the dimensionless time t, the rescaled 
Hubble parameter H is given by 

d 1 ' 2 1 /1.2 1-2 V /2 - - Y 2 Y 2 (T 2 

H = ^ = ^[& + ? +V <X>*)) ' nx,-) = (M 2 -^) 2 + ^ (12) 
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and M = 2.39 x 1(T 3 . 

We shall first examine the case where initially a ^> x, 1. Under these conditions, the last 
term in eq.([K]) dominates and this equation reduces to 

1 



X + 3H X + -a 2 x ^ . 



(13) 



Furthermore, the potential energy density V(x, &) i n ec i-(|l3) is dominated by its last term 
provided that yo ^ 2M 2 . Let us for the moment assume that, during the early stages of 
evolution, a remains almost constant. Then, for \y\ <C x<7/a/2, H is given approximately 
by 

1 „ „ 



H 



2v/3 



X° 



so that H and x satisfy the following differential equations 



x + ~y°xx + if x - o 



(14) 



(15) 



H + 3HH + -a 2 H ~ 
2 



Integrating once, with x(i = 0) = for simplicity, we get 



X(£) 2 ~x(0) 2 + -ln 



1 + — \ ( ' ) 
a 



' X(i) 



a 



(16) 



(17) 



H(ty ~H(oy + — In 



i + §-M 



-3*(f) 



(18) 



Using these equations one can show that x(i), x(t), and H(t) are initially (for %(£) > 0) 
monotonically decreasing functions of t. At some time to, becomes small enough and 
starts performing damped oscillations over the maximum at x — 0. For t >t Q , the continuity 
equation g = —3H(g + p) , with p = p/k 2 M^, p being the pressure, averaged over one 
oscillation of x becomes g = —?>H^g , where 7 = 1 for a x 2 potential @. This together with 
the fact that g is proportional to H 2 gives 
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H^i, t' = i-i + ^H(i )-\i>i ■ (19) 

The quantities x(t ) , H(to) as well as their time derivatives can be estimated as follows. 
When the x field starts oscillating, Hit) satisfies the equation 

H(t ) 2 ~ ~H{to) , (20) 

as one deduces from eq . (|T9|) . Eqs. fll4]) and (|i8|) then give 



Substituting into eq.(|20"D we, furthermore, get 



l-exp[-^(0) 2 ] 



(21) 



X(to) ^ ^#(*o) ^ ^^-^Ph^W 2 ] • (22) 

Note that eqs. ([21]) and (|22|) imply that |£| < for £ < i , as was required for the 

validity of eq. flT4]) . Also, for i < t , a ^> x ^ 1 an d in eq.(|TTJ) the Hubble parameter 
dominates over the "frequency" x/v^ . This equation then reduces to a ~ — x/\/3 and the 
variation of a till i is of order x(0)if(to) _1 ~ x(0)/<5" ^ 1- This justifies our assumption 
that a remains constant for i < to. 

Upon substituting eq.(|l^) into eq. (|i~3D we get 



For a constant, this is a differential equation of the Emden type whose solution satisfying 
the (approximate) boundary conditions at to is given by 



4 

X\t) ^ „, COS 



t > t • (24) 



Eq.(pT|), for t ^> toi averaged over one oscillation of x then gives 



t 6ot l 



This equation implies 



S 



where x(i) = t&(i), x(t) > O.We see that x(t) is a periodic function of time corresponding 
to a bounded motion of a; in a logarithmic potential. The turning point x m (> 0) is defined 
through the relation 

Unx(i ) 2 + x(t ) 2 = Unx 2 m (27) 
and eq.(p6|) becomes (x(t) > 0) 

_ _ ' (28) 



$L/ln^-)-$L/ln. 



. „ V37T 

where $(x) is the probability integral. The turning point x m can be estimated from eq. 
and the fact that x(i ) = a(i )i + a (to) ~ 0"(O). We get x m « x(t ) exp(3<r(0) 2 /8) > x(t ) 
for cr(0) ^> 1. This implies that the system, at t = t , starts from a position far away 
from the turning point which could subsequently be reached after the lapse of time Si ~ 
x mV / 3^/2v / 2 « (V3^/2V2) exp(3<r(0) 2 /8)x{i ) as one deduces from eq. (P^). Using x(to) = 
a(t )t « (2/3)a(0) J ff(t )" 1 « 2[1 - exp[-3x(0) 2 /4]]- 1/2 from eq.©, 5t can be related to 
the initial values of the fields x(0), cx(0). The time needed for the amplitude of x to drop to 
about 2M 2 /a, where the M 4 term starts dominating the potential energy V(x, &) in eq.(|T^), 
can be estimated from eq.(H) to be t D w 2/V3M 2 « 2 • 10 5 < 5t for <r(0) > 1 and £(0) not 
extremely small. We then deduce that a is an extremely slowly decreasing function of time, 
for t < t < t D , and, thus, our assumption that a remains constant is justified in this time 
interval too. 

For cosmic times > t D , the M 4 term dominates the potential energy in eq.([T2|) and the 
Hubble parameter becomes approximately constant and equal to H = M 2 /\/3 and remains 
so thereafter. Eq.(|T3"D still holds and, assuming that a remains essentially unchanged, \ 
performs damped oscillations of frequency a /y/2 ^> H . Eq. flTTl) averaged over one oscillation 
of x then becomes 



& + 3H&+-x 2 m a = , (29) 

where Xm is the amplitude of x- The decreasing "frequency" of a which is x m /2 < M 2 /& 
soon becomes much smaller than H and eq.(p9|) reduces to 

3ff*+^&7«0 • (30) 

The variation of a within one expansion time is then estimated to be \Aa/a\ ~ < 
M 4 /3(j 2 <C 1. This justifies the assumption that a does not change much. The overall 
conclusion is that, at times much larger than i D « 2/y/3M 2 , the x field falls into the valley 
of minima of the potential in eq. ( [12| ) and relaxes at the bottom of this valley whereas the a 
field still remains more or less unchanged. After that, the radiative corrections in eq.([5]) come 
into play and the system follows the valley of minima towards the super symmetric vacuum 
and, therefore, the hybrid inflationary scenario is realized for initial values of the fields 
satisfying the inequality a 1, x- These initial conditions, however, cannot be considered 
as totally satisfying because there is a considerable discrepancy between the initial values of 
the fields. Also, the inclusion of supergravity is expected to invalidate the above discussion 
of initial conditions which involve values of the field a > 1. 

For SUSY inflation based on radiative corrections to be considered as a fully successful 
inflationary scenario, one must show that it is obtained for a wide class of initial conditions 
which are more " natural" than the ones just discussed. This can be checked only numerically. 
To this end, we have conducted numerical integration of eqs.flTOD and ( pTTj ) for the following 
sets of initial conditions: 

1- 4>%>0.5, 6>6">0.5 with vanishing initial velocities. 

2- 1>X>0.1, 1><7>0.1 with vanishing initial velocities. 

3- l>x>0.1,l><r>0.1at the same points as in (2) introducing equal positive initial 
velocities so that the initial kinetic energy equals the corresponding potential energy. 
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The last two cases, although corresponding to pretty low values of the initial energy density, 
were included since our results there are expected to remain essentially unaffected if we 
replace global by local SUSY with canonical Kaehler potential. On the contrary, the results 
of case (1) should be strongly affected by supergravity corrections. To see this we reexamined 
this case by including supergravity corrections with minimal Kaehler potential. This is done 
by replacing the potential energy density in eq. (|i~2l) and its derivatives with respect to x an d 
a in eqs.fllO]) and fly]) by 



V(x, d)sa = exp 



2 



2 x + r 



+ (\x 2 - M 2 ) + ) (3D 



and its respective derivatives. 

The integration of the two coupled equations was performed by implementing a variant 
of the Bulrish-Stoer || variable step method in a Fortran program. As a general rule, the 
initial step for the dimensionless time variable t was chosen to be 2 while the sought accuracy 
was put to 1CT 15 . This choice was found to ensure reasonable stability in all cases. 

The results of our search are shown in Figs. 1,2,3, matching the corresponding sets of 
initial conditions, and 4 where supergravity corrections were included in the set of initial 
conditions (1). Each point on the a-\ plane corresponds to given initial conditions and the 
way it is depicted corresponds to a definite evolution type for the d-x system. We have used 
three different kinds of symbols. 

(a) Open circles: Both fields oscillate and fall towards the supersymmetric minima with- 
out producing any appreciable amount of inflation. 

(b) Filled triangles: Here, only x oscillates. The field a starts-off at large values (a > x) 
and remains almost constant for a very large period of time. The system eventually relaxes 
at the bottom of the valley of minima of V(x, <r) in eq . (|l"2"D . Its subsequent evolution along 
this valley produces an adequate amount of inflation. 
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(c) Filled circles: Both fields start oscillating at the beginning. Then a settles down at 
large values and the system follows an evolution of type (b) with adequate inflation. 
The evolution pattern (b) includes the limiting area a ^> 1, % already described and analyzed 
by means of semianalytic arguments. In case (c), although a starts mostly at moderate or 
small values a ~ Xi it appears to increase in amplitude absorbing energy from the fast os- 
cillating field x, creating eventually conditions favoring evolution of type (b). Here, we have 
an example of large energy transfer between two strongly-coupled non-linear oscillators. A 
closer look at Figs. 1-4 reveals that, in contrast to the case of the smooth hybrid inflation- 
ary model, connected regions corresponding to a definite evolution pattern do not appear to 
emerge except for the filled triangles area in Fig.l. This area, however, although leading to 
successful inflation cannot be considered as being "natural" since it requires relatively large 
differences between the initial values of the fields and is destroyed by supergravity as we 
shall soon see. Points following evolution type (c) are significantly better since, for many of 
them, the initial values of the fields are more or less of the same order of magnitude. These 
points, however, are scattered and do not appear to constitute an extended and well bounded 
continuous area as in the case of the smooth hybrid inflationary model. We have checked in 
more detail the area around most of the filled circle points and found out that almost all the 
points in their immediate neighborhood correspond to the evolution type (a) suggesting that 
evolution of type (c) occurs only for more or less isolated points. We conclude then that, in 
the model analyzed here, only about 10% of the examined points lead to successful inflation 
and lie in the area which can be considered "natural". Finally, the results with supergravity 
corrections included (Fig.4) show that the filled triangle area disappears and we are left with 
evolution of the types (a) and (b). This is to be expected, since the very steep walls of the 
potential in the a direction do not allow the a field to linger at large values for long. It is 
remarkable to realize that here also only about 10% of the examined points lead to adequate 
inflation. The initial energy density p can be calculated from eqs. fl9|), (|T2|) for global SUSY 
or (PD,(PT|) in the supergravity case with k given by eq.(|S|) where <tq/<7 c is taken close to 2 
@]. It turns out that, for many initial conditions corresponding to type (c) evolution, we 
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get adequately large initial energy densities. In summary, only a scattered set of " natural" 
but more or less isolated initial conditions, constituting about 10% of the examined cases, 
do lead to the SUSY inflationary model based on radiative corrections. It appears that this 
scenario is less satisfactory than smooth hybrid inflation which is obtainable for almost all 
initial conditions. 



ACKNOWLEDGMENTS 

This work was supported in part by E.U. grant ERBFMRXCT 960090 



13 



REFERENCES 

[1] A. D. Linde, Phys. Rev. D49 (1994) 748. 

[2] G. Lazarides and C. Panagiotakopoulos, Phys. Rev. D52 (1995) R559. 

[3] G. Dvali, Q. Shan and R. Schaefer, Phys. Rev.Lett. 73 (1994)1886 

[4] G. Lazarides, R.K. Schaefer and Q. Shan, Trieste preprint 10/96/115(1996) 
(Ihep-ph 96082561) 

[5] G. Lazarides, C. Panagiotakopoulos and N.D. Vlachos, Phys.Rev. D54 (1996) 1369. 

[6] For an earlier attempt at inflation with this superpotential, see E.J. Copeland, A.R. 
Liddle, D.H. Lyth, E.D. Stewart and D. Wands, Phys. Rev. D, 49, 6410 (1994). Su- 
persymmetric inflation has a long history and there are several reviews available, see, 
e.g., A.D. Linde, "Particle Physics and Inflationary Cosmology", Harwood Academic, 
Switzerland (1990). 

[7] M. Turner, Phys. Rev. D28 (1983) 1243. 

[8] W.H. Press, S.A. Teukolsky, W.T. Vetterling and B.P. Flannery, Numerical Recipes 
in Fortran: The Art of Scientific Computing, 2nd Ed. (Cambridge University Press, 
Cambridge, England, 1992) 



14 



FIGURE CAPTIONS 



Fig. 1: Evolution patterns for the G-\ system with vanishing initial velocities. Open circles represent points 
that do not lead to inflation. Filled triangles give adequate inflation having only the \ field oscillating. 
Filled circles give adequate inflation with both fields oscillating initially. 

Fig. 2: Same as in Fig.f . The initial values are now restricted to lie near the beginning of the axes. 
Fig. 3: Same as in Fig. 2. with positive initial velocities and equality of potential and kinetic energies. 
Fig. 4: Same as in Fig.f. Supcrgravity corrections are now being taken into account. 
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